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! Abstract 

> ■ 

00 ' One of the most dramatic discoveries made so far with the Rossi X-Ray 

Timing Explorer is that many accreting neutron stars with weak mag- 
netic fields generate strong, remarkably coherent, high-frequency X-ray 
brightness oscillations. The ~325-1200 Hz quasi-periodic oscillations 
(QPOs) observed in the accretion-powered emission are almost certainly 
produced by gas orbiting very close to the stellar surface and have fre- 
'■^ ' quencies related to the orbital frequencies of the gas. The ~360-600 Hz 

brightness oscillations seen during thermonuclear X-ray bursts are pro- 
duced by one or two hotter regions on the stellar surface and have 
\ frequencies equal to the stellar spin frequency or its first overtone. Mea- 

surements of these oscillations are providing tight upper bounds on the 
■ masses and radii of neutron stars, and important new constraints on the 

equation of state of neutron star matter. 
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1 Introduction 

Since the birth of X-ray astronomy 35 years ago, scientists have sought to use 
the X-radiation that comes from near the event horizons of black holes and 
the surfaces of neutron stars to probe quantitatively the strong gravitational 
fields near these objects and to determine the fundamental properties of dense 
matter (see, e.g., 0, ||). The Rossi X-Ray Timing Explorer (RXTE), which 
was launched on December 30, 1995, was specially designed to have the large 
area, microsecond time resolution, high telemetry bandwidth, and pointing 
flexibility needed to address these questions (see [pip . With RXTE, strong, 



high-frequency X-ray brightness oscillations have been discovered from at least 
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two black holes (see 0) and sixteen neutron stars (see [[13], p3[). As a result, 
we appear to be on the threshold of achieving this decades-old goal of X-ray 
astronomy. Here we focus on the oscillations discovered in neutron stars and 
their implications for the properties of these stars and for neutron star matter. 

2 High-Frequency Brightness Oscillations 

2.1 Observed Properties 

High-frequency X-ray brightness oscillations are observed both in the transient 
X-ray emission produced during type I (thermonuclear) X-ray bursts (see |[Lfc|| ) 
and in the persistent, accretion-powered X-ray emission (see fl23fl). 

The high-frequency oscillations observed during X-ray bursts have frequen- 
cies in the range 360-600 Hz and rms amplitudes at least as high as ~35% |19| . 



Only a single oscillation is observed during a burst, the oscillation appears to 
be highly coherent during the burst decay (see, e.g., [fTB]), and the frequency of 
the oscillations produced by a given star is always the same (measurements of 
the burst oscillations in 4U 1728—34 over about a year show that the timescale 
of any variation in the oscillation frequency is ^3000 yr [|17]]). The evidence 
is compelling that the oscillations seen during X-ray bursts are produced by 
regions of brighter X-ray emission that rotate with the star, and that the fre- 
quency of the burst oscillations is the stellar spin frequency or its first overtone 
(see |L3|, |19|). The spin frequencies of these neutron stars appear to be in the 



range 250-350 Hz [[T| . 



The kilohertz quasi-periodic oscillations (QPOs) observed in the persistent 
X-ray emission have frequencies in the range 325-1200 Hz, rms amplitudes as 



high as ~15%, and quality factors z//5z/ as high as ~200 | 2"5fl . Two kilohertz 
QPOs are commonly observed simultaneously in a given source (see Fig. 1). 
Although the frequencies of the two QPOs vary by hundreds of Hertz, the fre- 
quency separation Av between them appears to be nearly constant in almost 
all cases (see f23|). The frequency separation of the two kilohertz QPOs ob- 



served in a given star is closely equal to the spin frequency of the star inferred 
from its burst oscillations (see [|13|, |2( 

2.2 Origin of Kilohertz QPOs 



Although other types of models have been suggested [|, [22| , the fact that the 
separation Av between the frequencies of the two kilohertz QPOs observed 
from a given star is closely equal to the spin frequency of the star inferred from 
its burst oscillations is very strong evidence in favor of beat-frequency models 
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Fig. 1. Power spectrum 
of Sco X-l X-ray bright- 
ness variations, showing 
the two simultaneous kilo- 
hertz QPOs that are char- 
acteristic of the kilohertz 
QPO sources. The QPOs 
observed in Sco X-l are 
two of the weakest kilo- 
hertz QPOs detected, with 
rms amplitudes ~1%. The 
continuum power density 
is consistent with that ex- 
pected from photon count- 
ing noise alone. From 
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23]. In these models, the frequency of the higher-frequency QPO is the Kep- 
lerian frequency at a special radius and the frequency of the lower-frequency 
QPO is the difference between the higher frequency and the spin frequency 
^spin of the neutron star. The magnetospheric beat-frequency model, which 
was developed to explain the single, ~ 15-60 Hz X-ray brightness oscillations 
discovered in the Z sources a decade ago (see J5|), has been discussed as a 
possible explanation for the kilohertz QPOs. However, it is difficult to explain 
many basic features of the kilohertz QPOs using this model, including why 
there are two such QPOs (see [fL3|). 

The most fully developed and successful model of the kilohertz QPOs is 



the so-called sonic-point beat-frequency model This model is based on 

previous work ||, |K| which showed that the drag force produced by radiation 
from a central star can terminate a Keplerian disk flow near the star. In 
the sonic-point model, some accreting gas spirals inward in nearly circular 
orbits until it is close to the neutron star, where radiation forces or general 
relativistic effects cause a sudden increase in the inward radial velocity. The 
radius at which this occurs is conveniently referred to as the sonic radius, even 
though the transition to supersonic flow is not directly relevant in this model. 

The sharp increase in the inward velocity is usually caused by the drag 
exerted on the orbiting gas by radiation from the star, but may instead be 
caused by general relativistic corrections to Newtonian gravity, if the gas in 
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Fig. 2. View of the neutron star and accretion disk along the rotation axis of 
the disk, which is rotating counterclockwise in this view, (a) Spiral trajectory 
followed by a single element of gas as it falls supersonically from the sonic 
radius to the stellar surface, (b) Spiral pattern of higher gas density formed 
by gas streaming inward along spiral trajectories with the shape shown in 
(a), from a clump orbiting near the sonic radius. The trajectory and density 
pattern are from general relativistic gas dynamical calculations [13]. 



the Keplerian flow reaches the innermost stable circular orbit without being 
significantly affected by radiation. Gas streams inward from density fluctua- 
tions (clumps) orbiting near the sonic radius along tightly spiraling trajectories 
like that shown in Figure 2a, generating a more open spiral density pattern 
like that shown in Figure 2b. This pattern rotates around the star with the 
Keplerian orbital frequency at the sonic point, z/k s - Collision of the denser 
gas from the clumps with the stellar surface creates beams of brighter X-ray 
emission, like that indicated by the white dashed lines in Figure 2b. These 
beams move around the star's equator, generating a quasi-periodic brightness 
oscillation with frequency u^ s . Accreting gas is funneled to certain parts of the 
stellar surface by the star's weak magnetic field, producing weak X-ray beams 
that rotate with the star. These beams are overtaken by a given orbiting clump 
once each beat period, so the inward mass flux from the clumps, and hence 
the accretion luminosity, varies at the sonic-point beat frequency u Ks — f sp i n . 

The sonic-point model explains naturally why the separation between the 
frequencies of the two kilohertz QPOs is nearly constant and equal to the 
burst oscillation frequency or half this frequency. It is also consistent with 
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the accretion rates and stellar magnetic fields inferred previously and accounts 
for the main features of the kilohertz QPOs, including their high and variable 
frequencies, their high amplitudes and coherences, and the common occurrence 
of two simultaneous kilohertz QPOs Jl3 . 



3 Constraints from Kilohertz QPOs 

3.1 Nonrotating Stars 

In order to see how constraints on the equation of state of neutron star matter 
can be derived, suppose first that the star is not rotating and assume that, for 
the star in question, fQ P02 — the highest observed value of the frequency of the 
higher-frequency (Keplerian-frequency) QPO in the kilohertz QPO pair — is 
1220 Hz (this is the highest QPO frequency detected so far from any neutron 
star; see fll3l). Obviously, the orbital radius i? or b of the clumps producing the 
QPO must be greater than the stellar radius; R or ^ must also be greater than 
the radius R ms of the innermost stable circular orbit in order for the clumps to 
produce a wave train that lasts tens of oscillation periods, as observed. These 
requirements constrain the representative point of the star to lie in a pie-slice 
shaped region of the radius-mass plane (see Fig. 3a). This bounds the mass 
and radius of the star from above. In terms of ^^02 > these bounds are [K| 



M° ax = 2.2 (1000Hz/z/* PO2 ) M and < ax = 19.5 (1000 Hz/^ P02 ) km. (1) 

Figure 3b compares the regions of the radius-mass plane allowed for three 
values of ^qpo2 with the mass-radius relations for nonrotating stars given by 
five representative equations of state. 



3.2 Rotating Stars 

Rotation affects the structure of the star and the spacetime, altering the re- 
gion of the radius-mass plane allowed by a given value of ^qpo2- Slow rotation 
expands the allowed region whereas rapid rotation shrinks it greatly. The 
parameter that characterizes the importance of rotational effects is the dimen- 
sionless quantity j = cJ/GM 2 , where J and M are the angular momentum 
and gravitational mass of the star. For the spin frequencies ~300 Hz inferred 
for the neutron stars in the kilohertz QPO sources, j is ~ 0.1-0.3, depending 
on the mass of the star and the equation of state. For such small values of j, 
the structure of the star is almost unaffected and a treatment that is first-order 
in j is adequate. To this order, the existence of upper bounds on the mass and 
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Fig. 3. (a) Radius-mass plane, showing how to construct the region allowed for 
a nonrotating neutron star with Vqpq 2 = 1220 Hz. i? or b must be greater than 
the stellar radius, so the star's representative point must lie to the left of the 
(dashed) cubic curve M°(R orh ) that relates the star's mass to the radius of or- 
bits with frequency 1220 Hz. The high coherence of the oscillations constrains 
i?orb t° be greater than R ms , the radius of the innermost stable orbit, which 
means that the radius of the actual orbit must lie on the M°(i? or b) curve below 
its intersection with the (dotted) straight line M°(R ms ) that relates the star's 
mass to R ms . These requirements constrain the star's representative point to 
lie in the pie-slice shaped region enclosed by the solid line, (b) Comparison of 
the regions allowed for nonrotating stars with three different QPO frequencies 
with the mass-radius relations for nonrotating neutron stars given by five rep- 
resentative equations of state, (c) Regions allowed for rotating neutron stars 
with four values of j and Vqpq2 = 1220 Hz, when effects of the stellar spin are 
included to first-order (see text), (d) Illustrative Keplerian QPO frequency 
vs. accretion luminosity curve predicted by the sonic-point model. 
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Fig. 4. (a) Maximum frequency of a nearly circular orbit as a function of stellar 
spin frequency for equation of state FPS and stars with 1.4M (dotted line), 
1.6 Mq (dashed line), and 1.8 Mq (solid line) gravitational masses, (b) Maxi- 
mum frequency of a nearly circular orbit for a star of any mass, for the three 
equations of state indicated. See [13] for details and further discussion. 



radius can be proved analytically |y|. The bounds for prograde orbits are 



M max ^[l + 0.75jKin)]M° ax and R max « [1 + 0.20^*)]^ , (2) 

where Jospin) is the value of j for the observed stellar spin rate at the maximum 
allowed mass for the equation of state being considered and M^ ax and -R^ax 
are the bounds on the mass and radius for a nonrotating star (see eqs. [1]). 
Equations (2) show that the region allowed for a slowly rotating star is always 
larger than the region allowed for a nonrotating star, regardless of the equation 
of state. Figure lc illustrates the effects of slow stellar rotation on the allowed 
region of the radius-mass plane. Detailed calculations show that the mass 
of the neutron star in 4U 1636—536 must be less than 2.2 M & and its radius 
must be less than 17 km. The upper bounds may be smaller, depending on 



the equation of state assumed [[L3|, [14] 



For spin frequencies > 400 Hz , the structure of the star can be significantly 
affected, as well as the exterior spacetime. In this case both the structure of 
the star and the spacetime must be computed numerically for each assumed 
equation of state. Determination of the bounds on the mass and radius of a 
star for a given QPO frequency and equation of state therefore requires con- 
struction of a sequence of stellar models and spacetimes for different masses 
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using the equation of state, with i/ spin as measured at infinity held fixed |L3 



Such sequences have been constructed [[12| and show that if the neutron star 
is spinning rapidly, the constraints on the equation of state are tightened dra- 
matically. For example, if the ~580 Hz burst oscillation frequency observed 
in 4U 1636—536 is its spin frequency (rather than twice its spin frequency as 
indicated by the frequency separation of its two kilohertz QPOs), then very 
stiff equations of state like the tensor-interaction equation of state [15] are 
excluded by the 1220 Hz QPO already observed from this star. Regardless of 
the star's spin rate, a 1500 Hz QPO frequency would constrain the mass and 
radius of the neutron star to be less than ~1.7M and ~13 km, ruling out 



several equations of state that are currently astrophysically viable |L2[ . 

Figure 4a shows why slow rotation typically loosens the constraints on mass 
and radius implied by a given QPO frequency, whereas rapid rotation tightens 
them. At low spin rates the equatorial radius of the star is smaller than the 
radius of the innermost stable orbit and hence the maximum orbital frequency 
increases linearly with the star's spin rate (see eqs. [2]). In contrast, at high 
spin rates the equatorial radius becomes larger than the radius of the circular 
orbit that would be marginally stable, so no marginally stable orbit exists. 
The highest frequency orbit for a star of given mass is then the one just above 
the stellar surface, which increases in radius as the star spins faster, causing 
the highest possible orbital frequency to decrease with increasing spin rate. 
Figure 4b displays the maximum frequency of a nearly circular orbit for a star 
of any mass, for three different equations of state. Hence, if any neutron star 
is found to have spin and QPO frequencies that place its representative point 
above one of these curves, that equation of state is excluded. 



3.3 Innermost Stable Circular Orbit 

Establishing that an observed QPO frequency is the orbital frequency of the 
innermost stable circular orbit in an X-ray source would be an important step 
forward in our understanding of strong-field gravity and the properties of dense 
matter, because it would be the first confirmation of a prediction of general 
relativity in the strong-field regime and would also fix (for each equation of 
state) the mass of the neutron star involved. 

Given the fundamental significance of the detection of an innermost stable 
orbit, it is very important to establish what would constitute strong, rather 
than merely suggestive, evidence that the innermost stable orbit has been de- 
tected. Probably the most convincing signature would be a kilohertz QPO 
with a frequency that reproducibly increases steeply with increasing accretion 
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rate but then becomes constant and remains nearly constant as the accretion 
rate increases further. This behavior emerges naturally from general relativis- 
tic calculations of the gas dynamics and radiation transport in the sonic-point 
model (see Fig. 3d). The constant frequency should always be the same in a 
given source. Two other possible signatures of the innermost stable orbit are 



discussed in [13 



Several authors have recently suggested that innermost stable orbits have 
already been observed. Zhang et al. [27] suggested that the similarity of the 
highest QPO frequencies seen so far indicates that innermost stable orbits are 
being detected and that, based on the equations (1) for nonrotating stars, the 
neutron stars in all the kilohertz QPO sources therefore have masses close to 
2.0 M Q . Kaaret, Ford, & Chen M suggested that the ~800-900 Hz QPOs 



discovered in 4U 1608—52 |l| and 4U 1636—536 [26|, which were initially ob- 



served to have roughly constant frequencies, are generated by the beat of 
the spin frequency against the frequency of innermost stable orbits in these 
sources. However, no clear signature of an innermost stable circular orbit has 
so far been seen in any source. 

Indeed, subsequent observations of both 4U 1608—52 || and 4U 1636—536 



25] are inconsistent with the interpretation that the frequencies of the QPOs 



seen initially are related to the frequencies of innermost stable orbits around 



these stars. The 1171 Hz QPO seen by Zhang et al. [25], which was as- 
sumed by Kaaret et al. [[| to be at the frequency of the innermost stable 
orbit in 4U 1636—536 in order to estimate the mass of the star, was later seen 



at 1193 Hz [p5 and still later at 1220 Hz (W. Zhang, personal communica- 
tion). Hence, there is as yet no evidence for a maximum QPO frequency in 
4U 1636—536 and hence there is no basis for the suggestion that an innermost 
stable orbit has been seen in this source. A recent analysis of 4U 1608—52 
data by Mendez et al. || shows that this source has two kilohertz QPOs that 
vary with countrate just like the other sources. There is as yet no evidence for 
a maximum QPO frequency in 4U 1608—52 and hence there is no basis for the 
suggestion that an innermost stable orbit has been seen in this source, either. 



4 Constraints from Burst Oscillations 

As noted in §2, the strong (rms amplitudes up to at least 35%), high-frequency 
oscillations seen during X-ray bursts are thought to be caused by emission from 
a single or two nearly antipodal bright regions on the stellar surface, which 
produce large amplitude brightness oscillations at the stellar spin frequency or 
its first overtone as the star turns. The anisotropy at infinity of the radiation 
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Fig. 5. Upper bounds on the observed fractional rms amplitude of oscillations 
in the photon number flux during bursts, as a function of neutron star radius. 
The observer's line of sight is assumed to be in the rotational equator and the 
aberration and Doppler shifts caused by rotation are neglected, (a) Amplitude 
as a function of neutron star radius for isotropic emission (dotted line) or the 
peaked emission expected for a scattering atmosphere (solid line) from a single 
point, (b) Amplitude as a function of radius for the same intensity distribu- 
tions as in (a), but for emission from two identical, antipodal emitting points. 
See [jll| for details of the computational method and further discussion. 



emitted from such regions, and hence the amplitude of the burst oscillation, 
typically decreases with increasing gravitational light deflection by the star. 
Hence, in addition to any constraints on the radius of the star that may be 
derived from the X-ray spectra of the bursts, constraints on the compactness 
of the star (and hence the softness of neutron star matter) can be derived from 
the observed amplitudes of the burst oscillations [|TT], [lTj. Such constraints are 
particularly useful because they complement the constraints derived from the 
kilohertz QPOs, which constrain the stiffness of neutron star matter. 

Figure 5 shows results typical of general relativistic calculations of the 
maximum amplitudes of burst oscillations from one or two bright regions, as 
a function of stellar radius. Photon counting rate oscillations with 2-60 keV 
rms amplitudes as high as ~35% have been reported in 4U 1636—536 with 
a frequency equal to twice the apparent spin frequency of this neutron star. 
Figure 5b shows that an oscillation this large in the photon number flux would 
constrain the radius of this star to be greater than 5.5 M. However, the ampli- 
tude of the countrate oscillation measured at infinity can be larger or smaller 
than the amplitude of the oscillation in the photon flux at the star, depending 
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on the detector response as well as the angular dependence and spectrum of 
the emission from the stellar surface, the stellar spin rate, and whether there 
is any scattering material surrounding the star, so more detailed modeling will 



be required to extract the actual constraints on the compactness (see ['I 1 



5 Concluding Remarks 

The discovery using the Rossi X-Ray Timing Explorer that many neutron 
stars with weak magnetic fields produce strong ^300-1200 Hz X-ray bright- 
ness oscillations is a spectacular achievement that validates both the scientific 
expectations that led to the mission and the long years of hard work that were 
needed to bring it to fruition. The kilohertz QPOs discovered in the accretion- 
powered emission are already providing interesting new upper bounds on the 
masses and radii of neutron stars, and on the stiffness of neutron star matter. 
The high-frequency oscillations discovered in the emission during thermonu- 
clear X-ray bursts are likely to provide interesting new bounds on the com- 
pactness of neutron stars and hence on the softness of neutron star matter. 
Observation of a QPO with a frequency just 100 Hz higher than the highest 
frequency so far seen would exclude the stiff est proposed neutron star matter 
equations of state. 

Observation of innermost stable circular orbits would be the first confir- 
mation of a strong-field prediction of general relativity and would fix the mass 
of the star involved, for each equation of state considered. Although there is 
currently no strong evidence that an innermost stable circular orbit has been 
discovered around any of these neutron stars, there is reason to hope that such 
evidence may be forthcoming. Given the rapid pace of discoveries with RXTE, 
the prospects for obtaining compelling evidence of an innermost stable circular 
orbit appear good. 

This work was supported in part by NSF grant AST 96-18524, NASA 
grant NAG 5-2925, and NASA RXTE grants at the University of Illinois, and 
by NASA grant NAG 5-2868 at the University of Chicago. 



References 

[1] Berger, M., et al. 1996, ApJ, 469, L13 

[2] Epstein, R., Lamb, F.K., & Priedhorsky, W. 1986, Astrophysics of Time 
Variability in X-Ray and Gamma-Ray Sources, Los Alamos Science, 
No. 13 



11 



Kaaret, P., Ford, E. C, & Chen, K. 1997, ApJ, 480, L27 

Klein, R.I., Jernigan, J. G., Arons, J., Morgan, E. H., & Zhang, W. 1996, 
ApJ, 469, L119 

Lamb, F. K. 1991, in Neutron Stars: Theory and Observation, ed. J. 
Ventura & D. Pines (Dordrecht: Kluwer), 445 

Lamb, F. K., & Pines, D. 1979, Compact Galactic X-Ray Sources (Urbana: 
Univ. of Illinois Physics Dept.) 

McClintock, J. E. 1998, in Accretion Processes in Astrophysical Systems, 
ed. S. Holt & T. Kallman (AIP Conf. Proa), in press (|astro-ph/9802080|) 

Mendez, M., et al. 1998, ApJ, 494, L65 

Miller, M. C, & Lamb, F. K. 1993, ApJ, 413, L43 

-. 1996, ApJ, 470, 1033 

-. 1998, ApJ, in press (|astro-ph/9711325|) 



Miller, M. C, Lamb, F. K., & Cook, G. 1998, ApJ, submitted 

Miller, M. C, Lamb, F. K., & Psaltis, D. 1998, ApJ, in press ( |astro- 



ph/96091571 ) 
. 1998, ApJ, in preparation 

Pandharipande, V. R., & Smith, R. A. 1975, Nucl. Phys., A237, 507 

Smith, D.A., Morgan, E.,H., & Bradt, H. 1997, ApJ, 479, L137 

Strohmayer, T.E., talk presented at the 1997 November meeting of the 
High Energy Astrophysics Division of the American Astronomical Society 

Strohmayer, T. E., Swank, J. H., & Zhang, W. 1998, in The Active X-Ray 
Sky, eds. L. Scarsi, H. Bradt, P. Giommi, and F. Fiore, Nucl. Phys. B Proc. 
Suppl., in press flastr o-ph / 980 1 2 1 9 ) 



Strohmayer, T.E., Zhang, W., & Swank, J. H. 1997, ApJ, 487, L77 

Strohmayer, T., Zhang, W., Swank, J.H., Smale, A., Titarchuk, L., & 
Day, C. 1996, ApJ, 469, L9 

Swank, J., et al. 1995, in The Lives of Neutron Stars, ed. M. A. Alpar, U. 
Kiziloglu, & J. van Paradijs (Dordrecht: Kluwer), 525 

Titarchuk, L., & Muslimov, A. 1997, A&A, 323, L5 

van der Klis, M. 1998, in The Many Faces of Neutron Stars, Proc. NATO 
ASI, Lipari, Italy (Dordrecht: Kluwer), in press ( |astro-ph/9710016|) 



12 



[24] van der Klis, M., Wijnands, R., Home, K., & Chen, W. 1997, ApJ, 481, 
L97 

[25] Wijnands, R.A.D., et al. 1997, ApJ, 479, L141 

[26] Zhang, W., Lapidus, L, White, N. E., & Titarchuk, L. 1996, ApJ, 469, 
L17 

[27] Zhang, W., Strohmayer, T., & Swank, J. H. 1997, ApJ, 482, L167 



13 



